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We investigate the ground-state phase diagram of the Hubbard model for the ABjv_i chain 
with filling 1/N, where N is the number of atoms per unit cell. In the strong-coupling limit, a 
charge transition takes place from a band insulator (BI) to a correlated insulator (CI) for increasing 
on-site repulsion U and positive on-site energy difference A (energy at A sites lower than at B 
sites). In the weak-coupling limit, a bosonization analysis suggests that for N > 2 the physics 
is qualitatively similar to the case N — 2 which has already been studied: an intermediate phase 
emerges, which corresponds to a bond-ordered ferroelectric insulator (FI) with spontaneously broken 
inversion symmetry. We have determined the quantum phase diagram for the cases N = 3 and N = 4 
from the crossings of energy levels of appropriate excited states, which correspond to jumps in the 
charge and spin Berry phases, and from the change of sign of the localization parameter z c L . From 
these techniques we find that, quantitatively, the BI and FI phases are broader for N > 2 than 
when N — 2, in agreement with the bosonization analysis. Calculations of the Drude weight and z£ 
indicate that the system is insulating for all parameters, with the possible exception of the boundary 
between the BI and FI phases. 

PACS numbers: 71.10.Hf, 71.45.Lr, 78.20.Bh 



I. INTRODUCTION 

The half-filled Hubbard chain with alternating on-site 
energies ± tj A, known as the ionic Hubbard model (IHM), 
was proposed 1 ^ to describe the neutral-ionic transition 
in mixed-stack charge-transfer organic crystals such as 
tetrathiafulvalenc-p-chloranil.^4 4 > Interest in the model in- 
creased in the last decade due to its potential application 
to ferroelectric perovskitesi^SiL2i2iiSiiiii^ Although some 
details of the phase diagram, excitations, and expected 
physical properties of certain phases remain to be estab- 
lished, recent research has revealed the essential physics 
of the model, 10 ! 11 ! 12 . 13 ! 14 ! 15 ! 16 ! 17 ! 18 ! 19 

It is self-evident that in the strong-coupling limit (hop- 
ping t — > 0), the ground state is a band insulator (BI) 
for on-site repulsion U < A, when the sites with lower 
diagonal energy are doubly occupied, but is a type of 
Mott insulator (MI) for U > A, when all sites are singly 
occupied. However, for finite, small t, perturbational ap- 
proaches become invalid at U — A and non-trivial charge 
fluctuations persist even in the strongly coupled limit. 
Effective models around this limit have been proposed 
and analyzed 1 ^ but not yet studied numerically. Never- 
theless, following the initial proposal of Fabrizio et al.^ 
subsequent numerical studies f i2iiiii§ii2i and the obtain- 
ing of exact results for a closely related modelji 7 * 1 ^ it 
has become clear that the IHM chain has two transi- 
tions as U is increased. The first is a charge transition 
at U = U c , thought to be of the Ising typejiS^ii* 1 ^ from 
the BI to a bond-ordered, spontaneously dimerized, fer- 
roelectric insulator (FI). The second, when U is further 
increased, involves a vanishing of the spin gap at U s > U c 



at a Kosterlitz-Thouless transition between the FI and 
the MliiSaii The phase diagram has been constructed in 
full detail by following the crossing of appropriate ex- 
cited energy levels, which for this model turns out to be 
equivalent to the method of jumps in Berry phases (topo- 
logical transitions)^ For finite chains it has been shown 
that the topological transition at U c may be detectable 
in measurements of transport through annular molecules 
or nanodevicesiSi 

A particularly interesting feature of the model is the 
ferroelectric nature of the intermediate FI phase li 7 ^ 1 ^^ 
This phase results from an electronically induced Peierls 
instability, which generates a ferroelectric state with no 
ionic displacement. Interestingly, the elementary exci- 
tations of the FI phase have a fractional charge, which 
is proportional to the polarization^ Experimentally, a 
bond-ordered ferroelectric state has been observed in the 
pressure-temperature phase diagram of the prototypical 
compound tetrathiafulvalene-p-chloranili^ As noted by 
Egami et al.^ the microscopic origin of the displacive- 
type ferroelectric transition in covalent perovskite oxides 
such as BaTiOs remains unclear. 

The motivation for the present study is to investigate 
the changes in this type of phase diagram with the peri- 
odicity of the lattice, while retaining one electronegative 
ion and two electrons per unit cell, in such a way that 
for U = the system is clearly a BI. However, away from 
half-filling, i.e. for N > 2, the conventional Umklapp 
scattering term, which in the half-filled Hubbard chain is 
responsible for dynamical generation of a charge gap, is 
absent, and therefore one might expect metallic behav- 
ior or at minimum some qualitatively different physical 
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properties at larger values of U when N ^ 2. The model 
with N — 3 may also be relevant for the properties of 
doped, halogen-bridged binuclear metal chains (referred 
to as MMX chains) such as R 4 [Pt2(P205H2)4X.nH 2 0^ 
These compounds are related in turn to the quasi-one- 
dimensional MX complexes, which have been of interest 
for several decades in the broader context of the physics 
of electronic chain systems. 

The paper is organized as follows. In Sec. II we study 
the limit t — > to derive some results for the charge 
and spin gaps, and for the effective spin Hamiltonian. In 
Sec. Ill the opposite limit, of weak interactions, is treated 
by bosonization, and qualitative results are presented for 
the expected phases and correlation functions. Section 
IV contains a description of the numerical tools which 
are used in Sec. V to deduce the quantum phase diagram 
and to discuss the metallic or insulating character of the 
system. Section VI contains a summary and discussion. 

II. STRONG-COUPLING LIMIT 

We begin by expressing the Hamiltonian for the 
ABtv-i chain in a form to be used consistently through- 
out the following sections, 

H = -t ^2(c\ +la c la +H.c.) + U^n^nn 

+ ^A^ 4 . (1) 

i 

Here cj creates an electron at site i with spin a, nur = 
cJ^Cjo-, and rii = +n,-j i . The nearest-neighbor hopping 
amplitude is denoted by t, U is the on-site Hubbard in- 
teraction and A, = —A if i is a multiple of N and zero 
otherwise. Thus A > is the difference in on-site ener- 
gies between "metallic" (B or M) and "halogen" (A or 
X) sites on an ABat_i (Mn-iX) chain. 

A significant part of the physics of the ABat_i chain 
can be understood by considering the limit of small t. 
However, as in the case of the conventional IHM (N = 2), 
this limit presents particular complications for U — A. 
Although it is possible at this special point to make a 
canonical transformation retaining three states per site^ 
the resulting Hamiltonian is not trivial and no definitive 
conclusions may be drawn from it without further nu- 
merical analysis. This is beyond the scope of the current 
investigation, and in the remainder of this section we as- 
sume that \U — A| > t. 

For U < A, in a treatment discarding terms of order 
[i/(A — U)] 2 the A sites are doubly occupied and the 
B sites are empty. Thus the charges are ordered in a 
charge-density wave (CDW) as represented in Fig. ^a). 
The physics of this state can be described in terms of 
a one-particle picture in which the lower band is filled 
and the others are empty, yielding a BI. Again neglecting 
corrections of order t 2 , the charge and spin gap are both 
equal to A — U. We note that in contrast to the case 
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(a) 2k F CDW (LRO) 




(b)4k F CDW (LRO) + 2k p SDW 



FIG. 1: Charge and spin distribution in the strong-coupling 
limit, illustrated for an AB3 chain. 

N = 2, for N > 2 the charge distribution is altered 
radically if A is negative, because if t -C U the double 
occupancy at any site is no longer of order 1, but of 
order (t/U) 2 , but here we will not consider this situation 
further. 

Another simple limit is A = 0, where from the physics 
of the Hubbard model it is known that for any value of t 
the system is a MI for N = 2 but a Luttinger liquid for 
any N > 2&21& 

A further known result concerns the case N = 2 and 
U — A ^> t, where the system adopts a form of modified 
MI state, which we will call here the "correlated insula- 
tor" (CI) , with almost exactly one particle per site and a 
charge gap of approximately U — A. This state has gap- 
less spin excitations and long-ranged CDW order In 
spite of the fact that sites A and B are not equivalent, 
the system can be described by an effective Heisenberg 
model which is invariant under a one-site translation in- 
terchanging sites A and B&l^ However, the expectation 
values calculated with this effective Hamiltonian are not 
invariant under the same one-site translation due to the 
accompanying mapping of the operators^ The charge 
difference between A and B sites (the amplitude of the 
CDW) is at lowest order 22.2t 2 AU/{U 2 - A 2 ) 2 and the 
decay of the charge-charge correlation function with dis- 
tance d has the form d~ 3 ln~ 3 / 2 d for large d, with a pref- 
actor proportional to t 4 A 2 >i£*i& 

In order to extract analogous results for the case 
N > 2, we begin by considering the limit U — > +00, 
where the Hamiltonian for any N can be mapped to a 
non-interacting spinless system with an effective mag- 
netic fluxiSSiSi For N > 2 and any positive A, the sys- 
tem is an insulator with a charge distribution of the form 
depicted in Fig. db) . We will refer to this state as a CI 
to distinguish it from the charge distribution of the un- 
perturbed MI. The primary differences from the case of 
the CI for N = 2 are that a finite value of A is required 
to drive the system into an insulating state, and that the 
gap in this state is very small. For U — > +00 and N > 2, 
this gap is given by 

E c = 2i{cos[7r/A] - cos[2tt/A]} if A > t, 

E c = 2A/N if A< t. (2) 

For N > 2 and U - A > t but with finite U, the spin 
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degrees of freedom are important. A canonical transfor- 
mation which eliminates doubly occupied sites maps the 
model into a generalized t-J model, 



H = 



-t p(4 +la c ia + h.c.)p + A * n * 

ia i 

+ ti h P[c\ +s<J Ci-6^{2si-Si- 6 -\ni)]P 

iS=±l 

+ ^ Ji(si • s i+ i - |njn i+ i), (3) 



where P — Y[i(^~ n i't n il) ^ s ^ ne projector on the subspace 
of no double occupancy. The exchange parameter is given 
by 



Ji 



At 2 U 
U 2 - A 2 
At 2 



for i = lN or i = IN- 1, 



Ji = -jj- otherwise, 



(4) 



where / is an integer, and the correlated hopping term by 



±ch 



±ch 



ch 



u 
i 

2 V r 

t 2 



for i = IN, (I integer), 



U + A 



for i = IN + 1 or i = IN - 1, 



U 



otherwise. 



(5) 



As in the Hubbard model^ for the effective model of 
Eq. Q the charge degrees of freedom are independent of 
the spin of the particles in the limit U — > +oo, and can be 
described by a single Slater determinant corresponding to 
spinlcss fermions. The terms proportional to Ji and tf 1 
can be treated as a perturbation, and the resulting effec- 
tive Hamiltonian H s for the spin degrees of freedom takes 
the form of a Heisenberg model for a "squeezed" chain, 
by which is meant one in which the empty sites are elim- 
inated and the number of sites is equal to the number of 
electrons. The effective model H s is similar to general- 
ized t-J models which have been studied previously)22i31 
with one important exception: because Ji and tf 1 de- 
pend on the site, and the electrons which carry the spin 
are mobile (filling less than 1/2 for N > 2), the effec- 
tive exchange term depends on the charge configuration. 
However, because the effective hopping t is much larger 
than Ji and tf 1 for U — A ^> t (i.e. the charge velocity 
is much larger than the spin velocity), it is reasonable 
to assume that all spin degrees of freedom experience an 
average effective exchange interaction. 

Thus one may write 



Hs — ^ Jefi(Si 



(6) 



where 

Jeff 



N-l 



l^[J i (n i ni+i) / + 2tf(cJ_ 1 c i+1? 



i=0 



(7) 



and the expectation values (O) 1 of the operator O are 
evaluated in the spinlcss model. It is well known that 
H s has a gapless spectrum with algebraic decay of the 
spin-spin correlation function and only short-ranged an- 
tiferromagnetic order, as represented in Fig. ^b). 

In the effective Hamiltonian of Eq. © we have included 
only terms up to order t 2 which break the spin degeneracy 
of the U — > +oo limit. However, a next-nearest-neighbor 
exchange interaction J' appears at fourth order in t in the 
conventional IHM (N = 2), 2 and as U is decreased the 
quantity J 1 increases faster than does J c g. It is known 
that a spin gap opens for J' /J > 0.2411 ... in this spin 
modelj^i and thus a transition to a spin-gapped phase 
is expected as U is lowered^ One then expects that for 
N > 2 the presence of terms of higher order in t which are 
not included in Eq. © also drive the same sort of spin 
transition as U is lowered away from the strong-coupling 
limit. 

To summarize the results of this section, from a strong- 
coupling analysis of the ABat_! chain one may conclude 
that, for fixed, finite A and small U, in the limit t — ► 
the system has a BI phase. As U increases, the spin 
gap decreases within the BI regime and vanishes near 
U ~ A, as in the conventional IHM. For higher values 
of U the system adopts a CI phase which is similar to 
the MI in that charge degrees of freedom are high-lying 
and an effective spin model (albeit residing on a CDW 
background) describes the low-energy physics. However, 
from these considerations it is not possible to deduce the 
presence or absence of an intermediate FI phase. 



III. WEAK-COUPLING LIMIT 

In this section we analyze the model of Eq. using 
bosonization, a technique which is applicable in the limit 
(U, A) <C t. A conventional weak-coupling analysis im- 
plies a small interaction U, but in this case one may not 
expand around the non-interacting case (U — 0) with 
A > because in this regime the system is a BI and has 
no Fermi points. Instead, for U = A = the ABjv_! 
model reduces to a non-interacting chain with band fill- 
ing 1/N (two particles in N sites), for which the Fermi 
wave vector and Fermi velocity are given by 



Nao 



and vf = 2tao sm^aokp). 



(8) 



respectively, where do is the lattice spacing. 

As in the conventional IHM (N = 2), the bosonized 
expression for the A term [Eq. 124(1 ] is strongly rele- 
vant, and renders the usual renormalization-group treat- 
ment of the models invalid. For N = 2 some approx- 
imate and phenomcnological treatments have predicted 
the existence of the FI phase and its fractionally charged 
excitations>i2iii*i& In particular, by starting from the CI 
phase and integrating approximately over the charge de- 
grees of freedom, Fabrizio et al. predicted that the spin 
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transition takes place at a larger value of the on-site re- 
pulsion (U s ) than does the charge transition (U c ) i& Tak- 
ing the effective potential obtained from bosonization as 
a phcnomcnological Ginzburg-Landau free energy, these 
authors found that the excitations in the FI phase have 
a fractional charge which varies between 1 at the bound- 
ary with the BI phase and at the boundary with the 
CI phase. Thus the elementary excitation interpolates 
between an electron in the BI phase and a spinon in the 
CI phase, and its fractional charge is proportional to the 
electric polarization in the FI phase.' 1 — These excitations 
may be visualized as the topological excitations of an 
effective spin-1 chain for t -C U, AiiLi^ 

In the same spirit as the analysis of Fabrizio et ai, we 
identify the most important operators for N > 2 chains 
on the basis of their critical dimensions and their ex- 
pected effect on the physics. From these we infer the 
qualitative features of the phase diagram, which will be 
compared with numerical results in Sec. V. 

We linearize the spectrum and pass to the continuum 
limit by the substitution 



3 a 



^[e'^+^+e-^*^)] , (9) 



where the original lattice operators are decomposed into 
right- and left-moving components ^ + a (x) and a (x), 
and x = jag. Bosonization of these fields by a standard 
method for electrons with spin22i24 yields 



exp(=F# ±CT + 0±a), (10) 



where (j)+ a ((j>- a ) are right- (left)-moving Bose fields and 
9 r(y is the field dual to <j) ra . The Klein factors tv ra have the 
physical meaning of ladder operators which increase by 
one the number of fermions in branch r with spin a, and 
also serve to ensure that the anticommutation relations 
for electron fields of different spin are maintained. They 
are Hermitian and satisfy a Clifford algebra 3 ^ 



= {K-cr, K-a'} = 28 a 



(11) 



We define 



and introduce the linear combinations 



7#T 



V2 



ft 



(12) 



(13) 



to describe respectively the charge (<j> c ) and spin (cj) s ) 
degrees of freedom and their conjugate momenta. 

With the exception of the term in A^k F below, the 
above equations combined with appropriate operator 



product expansions give the bosonized expressions for the 
charge density 

p c (x) = :n iT +n a : -> -^(d x c/) c ) 

+A2k F cos(2kFX + V2(f) c ) cos V2<j) s 
-A 4k j,cos(4:k F x + V8(l)c), (14) 
and for the spin density 

Ps (x) = :|(n iT -n i; ): -> ^(d x (f> s ) 

+&2k F cos(2fci?x + y/2<j) c ) sin V2~4> s . (15) 

The term in A^k F , which is crucial in our analysis, is 
not present in a standard bosonization procedure but is 
generated in any model with density-density interaction 
terms. This may be seen by considering the lowest-order 
correction in the dressed charge operator, which includes 
terms of the form e~ ilkF "^\ ^ -,a, where each 
electron operator cj CT contributing to the four-fermion 
product introduces a factor of e~ lkp [see Eq. ©]. The 
presence of this term has been established definitively in 
the Hubbard model (any ABjv-i chain with A = 0). 36 

The coefficients ^ 2 (c F ! Aik F , and B2k F are non- 
universal parameters which depend on U and fulfil the 
conditions 



lim 

u^o 


A 2kF (U) 


= A = ^ 


lim 

£7^0 


B 2kF (U) 


- ±A9 

— 2-^1 ) 


lim 

u^o 


A 4kF (U) 


= o, 


lim 

U-> + oo 


A AkF (U) 


= Al 



(16) 



The coefficient A^k F has been calculated numerically for 
the particular case of the 0.55-filled Hubbard chaini 3 ^ It 
has a monotonic behavior, linear in U for U — > and 
with a downward curvature. For U — lOt it is of order 
0.12, already close to its saturation value A\ ~ 0.159 for 
U — > +oo. For the filling 1/N under consideration [two 
particles per unit cell, Eq. ©], a 2kp modulation cor- 
responds to the periodicity of the ABjv_i lattice, while 
Akp corresponds to a modulation of half this period in 
real space (Fig-QJ. 

To derive the bosonized expression of the on-site en- 
ergy term in the continuum limit, we note that it can be 
written in the form 



E A J n J = ~^ E E exp(2ik F a lj)r 

3 1=0 



N-l 



(17) 



By transforming the sum over j into an integral over 
x = jao, using Eq. (|14|l . and neglecting rapidly oscillating 
factors (which here means retaining only those terms in 
exp(ikpaon) with n a multiple of 2N), we obtain 



3 3 



A f 

-jjA 2 k F J dx cos(v / 20 c )cos(\/20 s ) 

A f 

+ jj-Aik F I dx cos(V80 c ). (18) 
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For the interaction term, with N > 2 one obtains the 
usual form 



nil; nil 



J dx | ^2 [{d x <j> c f - {d x <p s f 



U 



(27ra ) : 



cos 



(Vfys) . (19) 



Only for N = 2 does one have in addition the contribu- 
tion from Umklapp scattering processes 



B B 

(a) 



B B 

(b) 



FIG. 2: Bond-density distribution (a) in the BI phase and (b) 
in the FI phase of an AB2 chain. Thick lines correspond to 
high on-bond density, thin lines to less occupied bonds, and 
dashed lines to bonds with lowest occupation. 



u 



dx cos 



(20) 



which has the same form as the last term in Eq. I|18|l . We 
note that, even at weak coupling where the bosonization 
results are applicable, U AA/tk F and H\j m is the dom- 
inant term in the N = 2 caseiiS 

On including the non- interacting part, the final Hamil- 
tonian density for the bosonized model can be expressed 
as 



H e S = H c + H s + Has , 



(21) 



H c = 



H s 



■KVrKr 



M, 



2ttK, 



\d x ^ c f 



(27ra ) 

^nl(x) + ^-(d x ^ 



(2na ) 



■ cos 



(23) 



He 



na 



(v^ c ) cos (\/2> s ) . (24) 



Here Ti v is the moment conjugate to <p vi v c and v s are 
the charge and spin velocities, and K c is the charge- 
correlation exponent. For small interactions, the values 
of the different parameters are 

M c ~ U if N = 2, 

M c ~ A ■ Aik F (U) ifN>2, 

K c ~ 1 - Uao/nvF, v c ,v s ~vf 

M s ~ U, m cs ~A. (25) 

This Hamiltonian density coincides formally with the 
field theory studied previousljii2iii for the case N = 2 
corresponding to the IHM, with the sole difference aris- 
ing in the amplitude of the effective Umklapp scattering 
term M c . 

To discuss propeties of the different phases which ap- 
pear in the ABjv-i model, in addition to the on-site 
charge and spin density operators defined in Eqs. i|14|) 
and i|15|) we will use the bosonized expressions for the 



on-bond charge and spin density operatorS)2L2& 



^[H 2 c (x) + (d x <f> c f] 



+ cos((2i + 1)tt/7V + V20 c )cos V2</> s , (26) 



+ cos((2i + l) 7 r/A^ + \/20 c )sinv / 20 s . (27) 

To characterize the FI phase we introduce the order pa- 
rameter 

On = ^cos{2k F a j){B 3 - 



which represents the 2Uf Fourier component of the dif- 
ference between consecutive bond-density operators, and 
whose bosonized expression takes the form 



O-pi — sm{n/N) sm(V2ip c ) cos(V2(p s ). 



(28) 



We note that the operator O^i is antisymmetric with 
respect to inversion at the A sites (j — * —j). 

Classically, it is clear that minimization of the ionic 
term in Eq. I|24|l requires either (j> c = <j) s = or 
V2~4>c — V20S = 7T (mod 2tt). These values of <j> c and 
(j) s characterize the BI phase, as they ensure that bond- 
order parameter Opi and the site and bond spin densities 
[Eqs. (jT5|l and J27J] are suppressed, while there remains a 
long-ranged CDW order at 2k p [Eq. (|14f> ] with increased 
charge on the A atoms (for which x is multiple of Nao). 
The amplitude of the 4k p CDW is much smaller because 
it is proportional to Aik F - We stress that due to the 
symmetry of the Hamiltonian a 2kp CDW is present in 
all phases of the model. In contrast to the N = 2 case 
there is also a modulation of the bond order, or bond- 
order wave, in the BI phase with Bi oc cos[(2i + l)7r/iV] 
(|26[) as a consequence of the JV-site unit cell. This distri- 
bution of bond intensities is represented in Fig. 2(a), and 
is symmetric under reflection at the A sites (Ofi = 0). 

While \/8<fi c — (mod 27r) in the BI phase, minimizing 
the potential of H c [Eq. (|22p] requires y8(f) c = n (mod 
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2tt). Thus for fixed A there is a charge transition with 
increasing U for N = 2MU± For N > 2, M c is deter- 
mined not only by the on-site repulsion U but also by 
the ionicity parameter A. At small U, M c ~ AU/t, but 
for large U one expects M c ~ A as a result of the sat- 
uration of -44fc F (J7) in the Hubbard model (above). In 
consequence it is not possible to assure, as in the case 
N = 2, that this effect will overcome the ionicity for suf- 
ficiently large on-site repulsion, leading to a transition in 
the charge sector. However, from the results of the pre- 
vious section, in the strong-coupling limit (t — > 0) there 
is a charge transition for U = U c with U c ~ A. There 
is also a spin transition at U = U s because the BI phase 
has a spin gap which vanishes when U — A 3> t. Thus 
it appears that the same generic phase diagram indeed 
emerges for all values of N. 

For N = 2, the authors of Ref. ITol argued that U s > U c 
by describing the charge degrees of freedom in the CI 
phase in terms of a free boson of mass (charge gap) E c . 
On integrating over the charge degrees of freedom, the 
resulting effective Hamiltonian takes the form of Eq. ll'.'il) 
but with an effective renormalized mass 

M' s — M s ~ —8nt(A/E c ) 2 , (29) 

whence U s is determined by the condition M' s — 0. For 
U c < U < U s the system retains a spin gap as in the BI 
phase, with the values (p s — or tp a — tt/^/2 frozen, but 
the nature of the charge sector has changed. The spin 
correlations still decay exponentially due to the presence 
of the spin gap. 

An analysis of the intermediate regime may be per- 
formed for higher values of N by analogy with that 
applied to the model with N = 2. This type of 
treatmen1ji2*i& shows that the intermediate phase is 
characterized in the bosonized description by values of 
y84> c intermediate between and ir, and thus possesses 
fractional-charge excitations, broken inversion symmetry, 
and the non-zero polarizability of a ferroelectric phase. 
Because 4> c ^ 0, one observes from Eq. (|2"8l that the or- 
der parameter is finite, Ofi ^ 0, meaning that the charge 
distribution in the intermediate phase is characterized 
by broken inversion symmetry at the A sites. The long- 
ranged order of the BI phase, namely the 2fcp modulation 
of the site charge density and the TV-site-periodic modu- 
lation of the bond density, is retained. The distribution 
of bond charge density resulting from Eq. l|2t)|) for one 
of the two possible inequivalent choices of <j) v minimizing 
the energy is represented in Fig. 2(b), which makes clear 
the spontaneous breaking of parity in the FI phase. 

Qualitatively, the properties of the FI phase may be 
understood in this framework by taking the interac- 
tion part of the Hamiltonian density as a phenomeno- 
logical Ginzburg-Landau energy functional with effective 
interactions »i£ In the notation chosen here^ this energy 
functional takes the formic 

F = M c a 2 c + M s a 2 s + m cs a c a s , 
with a v — cos(v / 20„), (30) 



with M s < 0, M c > 2m cs , and all parameters duly renor- 
malized. The minima of F occur when a c = m C8 /(2M c ), 
a s = — 1, or a c = — to cs /(2M c ), a s = 1. A soliton be- 
tween two segments of the system characterized by these 
pairs of values corresponds to an elementary excitation 
of spin 1/2 and charge C = 1 — 2 arccos[m cs /(2M c )]/7r 
(i.e. proportional to the difference between the two val- 
ues of <fi c ). The polarization of the homogeneous system 
turns out to be ±eC/2ii& 

Finally, for U > U s , the spin gap is absent and the 
system is in the CI phase. With increasing on-site re- 
pulsion the amplitude of the 2kp modulations of the on- 
site charge density, A2k F , decreases and the <ikp compo- 
nent becomes dominant. Because of the spatial modula- 
tion of A there is long-ranged CDW order in the ground 
stat o 14 ' 15 as when N = 2. However, on top of this order 
the system also possesses the antiferromagnetic spin or- 
dering represented in Fig. ^b). As a consequence of the 
gapless character of the spin excitations, all correlation 
functions decay with a power-law form, including fluctu- 
ations of the charge density because of the strong charge- 
spin coupling [Eq. H24|)], These correlation functions have 
been considered by one of us for the case N = 

We conclude this section by extracting from the form 
of the bosonized effective Hamiltonian l|21l) the follow- 
ing differences between the well-characterized AB system 
and the ABat_i chain with N > 2. 

1) Because the amplitude of M c is very much smaller than 
for the AB chain, the BI gap decreases more slowly with 
increasing U and the BI phase should extend to larger 
values of U. 

2) The charge gap E c in the CI phase (U > U s ) is deter- 
mined not only by U but also by the ionicity parameter 
A. In particular, because the amplitude Aik F {U) of the 
Akp modulations saturates for large U, one expects that 
E c is determined solely by A and t for U — * +oo, in 
agreement with the results of the previous section. 

3) The smaller value of E c implies a larger renormaliza- 
tion of the effective spin mass M' s [Eq. I|29|) ]. and thus 
the FI phase is expected over a wider parameter range. 



IV. NUMERICAL METHODS 

We have performed numerical calculations of a num- 
ber of quantities characterizing the physical properties of 
ABjv-i chains for comparison with the analytical consid- 
erations of Sees. II and III, and to interpolate between 
the limits these represent. Here we describe the numeri- 
cal techniques applied to the Hamiltonian of Eq. . 

We have determined numerically the phase diagram for 
systems of L sites by considering the crossing of appro- 
priate excited energy levels which correspond to jumps 
in the charge (j c ) and spin (j s ) Berry phases, and by 
computing the charge (z£) and spin (z£) localization pa- 
rameters defined below. We have in addition calculated 
the Drude weight D c in some cases as a supplementary 
probe of metallic behavior. 
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The parameter j c is the Berry phase captured by the 
ground state for a ring threaded by a flux which is var- 
ied adiabatically from zero to two flux quantai^iSiil 
The spin Berry phase 7 S is the corresponding quantity 
for a situation in which oppositely directed fluxes are 
experienced by spin-up and spin-down particles42«ii4 
Specifically, under an applied flux (which we have scaled 
to be dimensionless) $> a hc/(2Tre) for spin a, the Hamil- 
tonian H [Eq. £[J] is transformed into a Hamiltonian 
H which differs from H in that the hopping term has 
the form -t E^cl+i^e 4 *^ + H.c). We denote by 

\g(<&-$,$l)) the ground state of i/($ T ,$J. The charge 
(spin) Berry phase 7 C (7^) is the overall phase captured 
by the state |ff($t, 3>j.)) on following adiabatically the cy- 



cle < $ < 2?r with $ 



$ ($ T 



Discretizing the interval < $ < 2n into N + 1 points 
$, r = 2nr/N (r — 0, N), the Berry phases are calculated 
from the numerically gauge-invariant expression 41 " 



N-2 



7c(s) = Jim Im{ln[TT(g($ r ,±$ r )| fl ($ r+ i,±* r+ i)) 

(31) 



r=0 



AT-.-OQ 

X< fl ($^_i,±* JV _i)|ff(27T,±27r))]}. 



The state |g(27r, 27r)) (1^(2^, — 27r))) is obtained from 
1 5 (0,0)) by applying the displacement operator £/£ (£/£), 



| ff (2^±2^))=[/f ) | 5 (0,0)), 



(32) 



with 



U C L = exp 



Uf 



exp 



2tt 



2tt „ 
i—D 
L 



, X = jao( n tf + riji), (33) 
3 

, D = J2ja a (n n - n n ). (34) 



Note that Eq. ll-'il't is simply a gauge transformation 
which transforms the ground state of H into that of 
H(2n, ±27r). The operator {/£, the exponential of the 
total-position operator, constitutes a translation in mo- 
mentum space which displaces all single-particle wave 
vectors by 2ir / Xi** Similarly, {/£ displaces the wave vec- 
tors of up and down particles in opposite directions. 

An important property of the charge Berry phase is 
that if the system is modified by some perturbation, the 
change in polarization Pf + is proportional to the cor- 
responding change in -f c Here P a is the contribution of 
electrons with spin a to the polarization of the system. 
Similarly, changes in j s are related to changes in the dif- 
ference Pj — P± between the electric polarizations for up 
and down spins ^ 



AP T ± AP l = eA7 c(s) /27r (mod e). 



(35) 



The ABjv-i chain has site inversion symmetry at the A 
sites for all N, and also at the B sites for N = 2. A 
crucial property for the purposes of the present analy- 
sis is that, in systems with site inversion symmetry, 7 C 



and 7 S can only take the values or tt (mod 2tt) (the 
argument of the logarithm in Eq. (|31|l becomes its own 
complex conjugate under inversion). Thus the Berry- 
phase vector 7 = (7c, 7 S ) cannot vary continuously, and 
a jump in 7 corresponds to a transition in at least one of 
the topological quantum numbers 7 c /tt and 7 s /7r. These 
topological transitions usually correspond to phase tran- 
sitions in the thermodynamic limit. As an example, in 
the strong-coupling limit the CI phase of the IHM has 
one charge at every site (1111. .. ), while the BI is char- 
acterized by a charge distribution of alternating electron 
pairs (2020. . . ). In order to transform from one to the 
other it is necessary to displace half of the charges by one 
lattice parameter, which corresponds to a change ±e/2 in 
the polarization of the system, and according to Eq. I|35|) 
the transition is accompanied by a jump of n in j c . Sim- 
ilarly, it can be shown that the opening of a spin gap in 
a Luttinger-liquid phase of a spin-rotationally invariant 
model is accompanied by a topological transition in 7^ 4£ 
Extrapolating the parameters for which these jumps oc- 
cur in systems with L < 16 to the thermodynamic limit 
has led to a very accurate determination of the quan- 
tum phase diagram of the Hubbard model with correlated 
hoppingiM44 Other successful applications of the method 
include the extended Hubbard chain with charge-dipole 
interactions 41 and the conventional IHM 12 (N = 2 of the 
present case). 

In the limit in which all charges are localized (generally 
the strong-coupling limit), the Berry phases are easy to 
calculate analytically. For t = one may choose a gauge 
in which all scalar products in Eq. (|3T|l are equal to 1, ex- 
cept possibly the last. The Berry phases are then defined 
by the argument of the exponential of in Eq. I|32|l . 
An illustrative example is the charge distribution of the 
BI for t — > [Fig. IHa)] in a system of M unit cells each 
of length N, such that L — MN: it is clear that D = 
(J33J and thus j s = 0, while [see Eq. 



X 



M-l 



NM(M-l), 



whence [/£ = exp[27ri(M — 1)] = 1 and 7 C = (mod 
27r). By contrast, for a system of particles with spin 
up localized at sites 0, N, 2N, . . . and those with spin 
down localized at positions N/2, N + N/2, 2N + N/2, . . . 
[Fig.Htb)], X = NM(M-l) + NM/2 and D = -ATM/2, 
leading to £/£ = U[ = — 1 and j c = j s = n; these values 
characterize the CI. 

The jumps in the Berry phases coincide in general 
with the crossing of pairs of low-lying excited states of 
a finite system. Thus the critical values U c and U s 
may be determined by the method of crossing excita- 
tion levels (MCEL), where this method is construed in 
a broader sense to be defined below. The relevant level 
crossing is found where the energy of the ground state 
|g($,±$)) as a function of the flux <!> reaches its max- 
imum valued This is usually equivalent to the situa- 
tion obtained by taking boundary conditions (BCs) oppo- 
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site to the "closed-shell" choice which leads to the mini- 
mum energy^ 2 " Identifying the two crossing levels by their 
quantum numbers and finding the parameters for which 
their energies coincide is naturally equivalent to locat- 
ing the jump in the Berry phase, and the former proce- 
dure saves computer time by avoiding the evaluation of 
Eq. I(31|) . Thus we have followed this method in deter- 
mining the majority of the data to be shown in Sec. V, 
although we have also used Eq. (|31() for verification of 
our results. 

In its more restricted sense, the MCEL is based on 
identifying the appropriate levels by conformal field the- 
ory with renormalization-group analysis This is 
a weak-coupling approach which takes advantage of the 
fact that in conformally invariant systems of finite size 
the smallest excitation gap corresponds to the domi- 
nant correlations at large distances. In previous stud- 
ies, which included the Hubbard model with correlated 
hopping^44i2& the extended Hubbard modelfi24& and 
the conventional IHM, 12 the restricted MCEL and the 
jumps of Berry phases were found to coincide, giving sup- 
port to the method from both weak- and strong-coupling 
analyses. 

From the arguments presented in the strong-coupling 
limit, it is evident that the jump in 7 C characterizes a 
sharp transition with a reordering of charge. One might 
expect by continuity that this jump characterizes the 
charge transition also at weak coupling. While in previ- 
ous studies (above) the jump in j c does coincide with the 
crossing of appropriate excitations determined from field- 
theoretical arguments for larger N it is not known 
whether the level crossing which corresponds to the jump 
in 7 C has support from a weak-coupling approach. In the 
following, the MCEL refers more broadly to a crossing 
which corresponds to the jump in a Berry phase, but in 
the case of j c this correspondence is not justified by field- 
theoretical methods. When considering 7 S in any model 
with SU(2) symmetry, this crossing always coincides with 
that corresponding to the Kosterlitz-Thouless transition 
for the opening of a spin gap according to conformal field 
theoryia- 

The specific crossing levels in the MCEL are deter- 
mined by the quantum numbers specifying their parity, 
spin, and total momentum, plus the BCs for which the 
crossing occurs. The spin transition is determined from 
the crossing of an even singlet and an odd triplet excited 
state, both with total momentum K = 0, when working 
with open-shell BCs. This is the crossing which corre- 
sponds to the jump in the spin Berry phased and has 
support from field-theoretic arguments £LML The charge 
transition is determined from the crossing of an even and 
an odd singlet, both with K — n/ao, again calculated 
with open-shell BCs. This crossing signals the jump in 
the charge Berry phase, which as explained above denotes 
a reordering of charge leading to a jump in the polariza- 
tion of the system. However, as in other recent calcula- 
tions in which crossing of levels was used to determine 
phase transitions in two-dimensional spin systems^*^ 



a systematic demonstration that this procedure should 
remain valid in the weak-coupling limit is lacking. For 
this reason, and because of other shortcomings described 
below, we reinforce the results obtained from the MCEL 
by comparison with those obtained from the localization 
operators. 

The localization operators are obtained from the ex- 
pectation values of the displacement operators 

z «») = {g \U<>) ]g) . (36) 

\z1\ was first proposed by Resta and Sorella as an in- 
dicator of localization in extended systems^ By appeal- 
ing to symmetry properties, Ortiz and one of the au- 
thors demonstrated that for translationally invariant in- 
teracting systems with a rational number n/l of parti- 
cles per unit cell the correct definition of the matrix ele- 
ment is (<?|(?7£) |<?)i" This definition was used to charac- 
terize metal-insulator and metal-superconducting tran- 
sitions in one-dimensional lattice models4^£*£S In the 
thermodynamic limit, L — > oo, |z£| is equal to 1 for 
a periodic metallic system and to for the insulating 
state. This result holds also for non-interacting disor- 
dered systems^ For an intuitive understanding we note 
that in a metallic system with a well-defined Fermi sur- 
face, £/£ displaces the Fermi surface by a wave vector 
2ir/L, and thus (g\Ul\g) = 0. By contrast, in a non- 
interacting band insulator C/£|<?) coincides with \g) and 

4 = 1- 

The determination of #2 in the thermodynamic 
limit provides a means of calculating the Berry 
phasea 8 i 42 i 44 i 45 . 50 

7c(a) = r lim Imlnz2 (s) . (37) 

For finite L, the right-hand side of Eq. I|37|) is actually 
equivalent to Eq. i|31|) if only a single point is used in 
the discretization of the flux. Thus it is not surpris- 
ing that direct calculation of the Berry phase (or the 
MCEL) leads to a more precise extrapolated result than 
does Eq. l(3T|) .^Jk For comparison with these two types 

of analysis we have performed calculations of z c L , as 
it emerges that these are important in determining the 
phase diagram of the AB3 model at low values of A. In 
addition, the displacement operators have a direct rela- 
tionship with the bosonic fields introduced in Sec. Ill: it 
was shown recently that J7£ and U[ can be constructed 
as exponentials of the average charge and spin fields ^ 

Ul = cxp(zV8C), € = \ J dxcj> v {x). (38) 

In the region of parameters such that the interaction M c 
in Eq. 1)22(1 is much greater than m C8 [Eq. 1(24(1 ] . (f> c (x) be- 
comes locked at the value which minimizes M c cos(v / 80 c ), 
leading to V8<p^ = j c — n (mod 2n). If instead m cs 
is dominant, s/%<j)% — 7 C = 0. Similarly, by mini- 
mizing the potential for large positive M' s one expects 
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FIG. 3: Critical value of on-site repulsion at the charge transi- 
tion as a function of 1/L 2 , determined by the MCEL (circles) 
and by the condition z£ = (squares) for the AB2 chain with 
A = 2. 
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FIG. 4: Critical value of on-site repulsion at the spin transi- 
tion as a function of 1/L 2 , determined by the MCEL for the 
AB2 chain with two values of A. 



= 7„ = 7r (in fact M' s renormalizes to zero when 
it is positive, but the result 7 S = tt is robust), while 
for negative M' s one has instead V8(f>s — 7 S = 0. This 
is in agreement with the values obtained in the strong- 
coupling limit for the BI and the CI. 

In the intermediate FI phase the spin gap remains 
open, and thus the spin Berry phase 7^ is zero, as in the 
BI. The situation with the charge Berry phase is more 
delicate: in the thermodynamic limit the inversion sym- 
metry is broken spontaneously, leading to intermediate 
values of the polarization and also to fractional values of 
V8<Pc, as has been shown in detail for N = 2^ Both 
quantities are also related to the charge of the fractional 
excitations of the modeliiS* 1 ^ However, for all finite sys- 
tems, and in particular those of our numerical analysis, 
the ground state of the FI retains inversion symmetry 
and 7 C = ir as in the CI phase. 

In summary, as for the case N = 2^ the Berry phase 
vector 7 has the value (0, 0) in the BI phase, (ir, 0) in the 
FI phase and (it, tt) in the CI phase. 

Finally, to investigate the possibility of a metallic 
phase near the charge transition, we have also calculated 
the Drude weight 



D r = - 



Ld 2 E(L,<S>) 
2 



<9$2 



l*=*o 



(39) 



where E(L, $) is the energy of the ground state for a sys- 
tem of length L with (anti-)periodic BCs in the presence 
of a flux &hc/ (2-7re) threading the ring. For the calcula- 
tion of Ul {s) and D c we have used BCs (or equivalently 
the choice of $0) corresponding to the "closed-shell" sit- 
uation which gives the minimum in the energy as a func- 
tion of flux. 



V. RESULTS 

To determine numerically the phase diagram of the 
AB 2 (N = 3) and AB 3 chains (N = 4) we have used both 
the MCEL, which corresponds to the jumps in charge 
and spin Berry phases (above), and the sign of z£, both 
calculated by exact diagonalization with the Lanczos al- 
gorithm. The largest accessible system sizes are then 
L = 15 for N — 3 and L = 16 for A" = 4, which are 
sufficient to allow a finite-size scaling analysis of limited 
accuracy. The two calculational methods should give the 
same result in the thermodynamic limit [Eq. Q37[l]. al- 
though as discussed in the previous section the former 
is usually more accurate^ However, we have found that 
for A < t the level crossings have large finite-size effects, 
particularly for A" = 4, which limits the comparison of 
the two techniques. We have extrapolated the transi- 
tion parameters using a quadratic polynomial in 1/L 2 , 
a dependence expected from conformal invariance and 
generally used in the MCEL. 31 i 46 i 47 

Examples of the size-dependence of the resulting criti- 
cal values at the charge (U c ) and spin (U s ) transitions are 
shown in Figs. and ^respectively. The unit of energy is 
chosen as t = 1. For U c we have also tested an extrapo- 
lation with 1/L 3 dependence^ but the results were sig- 
nificantly inferior and displayed larger errors. For A > 2 
the error in the extrapolated value of U c obtained from 
the condition z£ = is an order of magnitude larger than 
that obtained from the MCEL. Thus the MCEL may be 
taken to provide the more accurate results for A > 2, 
as expected. For the majority of the parameter range 
these methods extrapolate acceptably well, and in a sim- 
ilar manner, towards a fixed value. However, for smaller 
values of A, i.e. in the weak-coupling regime, the quality 
of the fit deteriorates rapidly for both methods, but par- 
ticularly for the MCEL. As discussed below, this method 
becomes unreliable for N — 4 and A < 1. 

Figure0]shows the value of the spin transition (U s ) in 
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FIG. 5: Scaling dimensions as a function of [7 for the AB3 
chain with L = 12 and A = 2. 
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FIG. 6: Phase diagram of the AB2 chain. 



the AB2 chain for different system sizes, obtained using 
the MCEL. Again the data extrapolate in a satisfactory 
manner for larger values of A, but a comparison of the 
two curves shows that the quality of the 1/L 2 fit for the 
size dependence of U s also deteriorates with decreasing 
A. A test of the MCEL based on conformal invariance 
and the renormalization group, and of the Kosterlitz- 
Thouless character of the spin transition, can be made 
by calculating the scaling dimensions x a ,\ of the singlet 
and of the triplet excitationsi^iilii According to the 
predictions of conformal field theory, the excitation en- 
ergies should scale according to 

A^ CTii S -j^x^i, (40) 

where A.E a< i (AE at 2) is the difference between the energy 
of the lowest singlet (triplet) state with BCs opposite 
to the closed-shell situation and the ground-state energy 
calculated with closed-shell BCs. The spin velocity v s is 
calculated from 

v, - S^^jM , ,41) 

where E(S, K) is the lowest energy in the sector of total 
spin S and total wave vector K. 

The values of and x a ^ cross at U s . Their aver- 
age value x av — (x<r,i + 3x ff ,2)/4 should be equal to 1/2 
near the crossing, and for U > U s where the spin gap is 
closed. Numerical results for the scaling dimensions are 
shown in Fig. [SJ Clearly x av is indeed close to 1/2, sup- 
porting the validity of the MCEL technique for the spin 
transition. For the reasons discussed in Sec. IV a similar 
analysis at U c is less reliable. Systems of larger size may 
be studied only by the density-matrix renormalization- 
group technique, which would offer a means of refining 
the quantitative aspects of the results below, as well as 
the possibility of investigating the critical behavior at the 
two transitions. 



Having demonstrated the application of the numerical 
methods discussed in Sec. IV, we are now in a position to 
discuss the phase diagrams of ABjv-i chains. In Fig. El 
we show the extrapolated values of U c — A and U s — A 
as functions of A for the AB2 system. The chain lengths 
used were L = 6, 9, 12, and 15. For U c we compare 
the results of the MCEL with those obtained from the 
change of sign of z£ [cf. Fig. |J3}] . The difference between 
the two methods is significant, particularly for A < 1, 
where both methods have larger errors. For A < 0.75, 
the points displayed for U s , and for U c obtained from z£, 
were calculated by fixing U and extrapolating the critical 
value of A. In comparison with the results obtained by 
fixing A, this procedure shifts the curves to smaller values 
of A, particularly for the MCEL. For the reasons stated 
above, at least for A > 2, we take the MCEL results to 
be more accurate than those obtained from the change of 
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FIG. 7: Drude weight and z c L as functions of U in the AB2 
chain with A = 10 for a range of system sizes. The criti- 
cal interaction strength obtained from the MCEL and by ex- 
traploating to the thermodynamic limit is U c = 12.94 for this 
value of A. 
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FIG. 8: Phase diagram of the AB3 chain. 



sign of z£. However, it is clear from Fig. that the two 
methods agree on a semiquantitative level. 

As expected from the strong-coupling analysis of 
Sec. II, the differences U c — A and U s — A are of the or- 
der of the hopping t. For comparison, we have drawn two 
horizontal lines which correspond to the charge and spin 
transitions obtained with the MCEL for the AB chain 
(the conventional IHMi£) for t -C U, A: U c - A = 1.334 
and J7 S — A = 1.914. For the AB2 chain J7 C is shifted to 
larger values, and also the difference U s — U c becomes 
larger for A > t. Thus while there are no qualita- 
tive changes to the phase diagram, the widths of the BI 
regime and also of the intermediate FI phase increase on 
passing from AB to AB2. This is in agreement with the 
conclusions of the bosonization treatment presented in 
Sec. III. For A < t, the numerical results show a ten- 
dency towards a negative difference U s — U c ; this result 
is difficult to justify on physical grounds, such as those 
underlying the bosonization treatment, and is in all prob- 
ability an artifact of finite-size effects, which are of the 
order of \U S — U c \ when U s — U c is negative. 

Figure0shows the dependence of the Drude weight D c 
and of the expectation value of the localization operator 
z£ as functions of U for A = 104 and for all system sizes 
studied. The evolution of these quantities with L (D c 
decreasing and |z£| increasing) for U 7^ U c indicates that 
the system is insulating over the entire range of param- 
eters. However, because the system sizes are small one 
may not exclude metallic behavior at or near the point 
U = U c . The vanishing value of |z£| at this point and 
the peak in D c are consistent with this possibility, but 
also simply with a larger localization length for U ~ U c . 
From our theoretical and numerical analyses it is not pos- 
sible to establish whether or not D c extrapolates to a fi- 
nite value, which would indicate metallic behavior in the 
thermodynamic limit. 

In Fig. [8] we show the extrapolated values of U c — A 
and U s — A as functions of A for the AB3 chain. The 
system sizes used for this case were L = 8, 12, and 16, 
combined with a limited number of runs with L = 20. 
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FIG. 9: Drude weight and z% as functions of U in the AB3 
chain with A = 10 for a range of system sizes. The criti- 
cal interaction strength obtained from the MCEL and by ex- 
traploating to the thermodynamic limit is U c = 13.20 for this 
value of A. 



For A < 1, all of the points displayed except those for U c 
obtained within the MCEL were calculated by fixing U 
and extrapolating the critical value of A. Surprisingly, 
U c obtained from the MCEL (or the jump in the charge 
Berry phase) appears to diverge as A — > 0. To be specific, 
for U — 50 (not shown) the extrapolated value of A at 
the transition is A c ~ 0.6, while for L = 8, A c ~ 3; these 
results reflect the presence of very large finite-size effects. 

This type of behavior is completely inconsistent with 
the strong-coupling results presented in Sec. II, which 
show that for any A > and U — A ^> t, the system is 
a CI with a charge distribution of the form represented 
schematically in Fig. ^b). By contrast, the transition 
obtained from z° h is consistent with the strong-coupling 
results. This fact, taken in combination with the very 
large finite-size effects of the MCEL for small A, sup- 
port the conclusion that the level crossing which is fol- 
lowed fails to describe the charge transition for A < 1. 
This difficulty may be related to the loss of translational 
symmetry in comparison with the Hubbard model. In 
fact, in the example of the extended Hubbard model 
with correlated hopping, where a similar charge transi- 
tion occurs^i2iiii& the quantum numbers for the ap- 
propriate crossing energy levels can be identified from 
a weak-coupling analysis and correspond to wave vector 
K = 7r/ao4& When the one-site translational symme- 
try is lost, this wave vector is mixed with others and the 
Lanczos method, which captures the lowest eigenstate for 
the corresponding wave vector in the reduced Brillouin 
zone, no longer follows necessarily the correct excitation 
in the weakly coupled case. This problem arises also in 
the extended IHMf±2 

We proceed by disregarding the points obtained for the 
charge transition with the MCEL for A < 1, and in this 
regime adopt the values obtained from z£. The phase 
diagram of the AB 3 chain (Fig. [SJ is very similar to that 
of AB2. Only at the quantitative level is there a moder- 
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ate shift of U c and U s to still larger values and a small 
increase in the width of the region occupied by the FI 
phase. In Fig.|Hlwe show also the values of U s for the spin 
transition obtained from the condition z s L = 0. These 
are in very close coincidence with the results obtained 
from the MCEL for A > 2; for small A both methods 
are again plagued by large finite-size effects, but the dis- 
crepancy between the extrapolated values is significantly 
smaller than in the case of the charge transition. 

Finally, Fig.[5]displays the dependence of D c and z c L on 
U in an AB 3 chain with A = lOi for several system sizes. 
The results are qualitatively and quantitatively very sim- 
ilar to those obtained for the AB2 system (Fig. [7J . At 
U > U c some indications of a smaller gap and longer lo- 
calization length might have been expected as the size 
of the unit cell N increases, as suggested by the results 
of Sec. II, but in this respect no significant differences 
are found between the AB2 and AB3 chains investigated 
here. 



VI. SUMMARY AND DISCUSSION 

We have performed a systematic analysis of the gen- 
eralized ionic Hubbard model (IHM) in one dimension. 
Motivated by the complex behavior emerging in the AB 
chain, where the commensurate filling ensures a domi- 
nant role for Umklapp scattering terms, we have inves- 
tigated the phase diagrams for ABjv-i chains with unit 
cell size N and filling 1/N. We have employed analytical 
considerations in both weak- and strong-coupling limits, 
and sophisticated numerical techniques based on Lanczos 
exact diagonalization, applying these specifically to the 
cases N = 3 and N = 4. 

Despite the absence of direct Umklapp terms away 
from N = 2, we find that the qualitative features of the 
AB7V-1 phase diagram are essentially identical to those of 
the conventional IHM. The generalized model always pos- 
sesses three phases as a function increasing values of U: 
a band insulator (BI) for small U/A, a correlated insula- 
tor (CI) for large U/A which is a form of Mott insulator, 
and between these an intermediate ferroelectric insula- 
tor (FI) with spontaneous breaking of inversion symme- 
try and fractional excitations. The generic nature of the 
phase diagram is explained within a bosonization anal- 
ysis by the emergence of higher-order interaction terms 
related to the charge-density modulation caused by the 
ionic term. 

The general similarity of all the models in this class 
masks some important differences. Qualitatively, the CI 
phase for large U is not strictly a Mott insulator because 
the number of particles is less than the number of sites. 
As a quantitative consequence, the localization length is 
therefore larger than one lattice parameter and the gap is 
much smaller than in the conventional (N = 2) IHM. In 
the latter, for large U the gap is of order U — A, whereas 
for N > 2 it is given by Eq. (J2J) and vanishes for A = 0. 
More generally, the relatively weaker Umklapp terms re- 



sult in larger values of U being required to compete with 
the ionic term A, which favors the BI phase, and thus the 
FI and CI phase boundaries are pushed to larger U with 
increasing N; both BI and FI regimes are broadened. 

We have implemented the numerical calculation of a 
number of characteristic quantities which provide valu- 
able insight into the properties of the ground and ex- 
cited states. For the two phase boundaries we have em- 
ployed the method of crossing excitation levels (MCEL), 
which is equivalent to following discontinuous steps in 
the charge and spin Berry phases, quantities which sig- 
nal the presence of a transition even on a small system. 
This calculation is compared with the results obtained 
from the vanishing of the expectation value of the dis- 
placement operators, which is found to be a less accu- 
rate but sometimes more reliable technique. We have 
also computed the Drude weight with a view to testing 
the possible metallic nature of the system at the phase 
boundaries between the insulating phases. 

While the numerical analysis serves as a useful test of 
these techniques, and allows us to make certain power- 
ful qualitative statements, the results are also subject to 
strong finite-size effects. In particular, the MCEL has 
significantly larger finite-size effects for N > 2 than at 
N = 2, due probably to the much smaller gap in the CI 
phase and the relatively lower translational symmetry. 
Indeed, at the charge transition for A < t the phase of 
z c L gives more reliable information. However, even with 
the restricted system sizes accessible in our study, the 
majority of the results obtained in Sec. V can be said 
to be of semi-quantitative accuracy. We comment that 
further studies on systems of larger size are possible us- 
ing the density-matrix renormalization-group technique: 
with the aid of the quantities we have defined and an- 
alyzed in Sees. IV and V, such calculations might con- 
firm our findings with greater accuracy, establish more 
precisely the boundary of the charge transition, and pos- 
sibly also permit the analysis of critical behavior at the 
two transitions. 

In the same way that the AB chain has for many years 
been of interest in the context of quasi-one-dimensional 
organic MX complexes, our considerations concerning 
the AB2 chain are expected to be relevant for the prop- 
erties of MMX chains. These are halogen-bridged bin- 
uclear metal chains such as RafPta^aOsHa^X.nHaO, 
on which experimental studies would be welcomed. We 
comment that the greater extent in parameter space of 
the FI phase for N > 2 systems may even be important 
for technological applications. Finally, with regard to the 
possible metallic properties of ABjv-i systems, we have 
found that despite being at fillings 1/N away from the 
half-filled situation, these remain insulating for all pa- 
rameter regimes other than A = 0, as a consequence of 
the commensurate nature of their electron number. Only 
at the charge transition do we obtain evidence from the 
Drude weight and from the charge localization operator 
of a possible metallic point in the phase diagram. 
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